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Abstract. A system of two coupled nonlinear Schrodinger equations is treated. In ad- 
dition, a linear coupling which models an external driven field described by the Rabi 
frequency is considered. Asymptotics for large Rabi frequency are carried out. Conver- 
gence in the appropriate Strichartz space is proven. As a consequence, the global existence 
of the limiting system gives a criterion for the long time behavior of the original system. 



1. Introduction 
We consider the following system of nonlinear Schrodinger equations 



1. , , i\ 

-\x 

idt1p2 = -^A?/;2 + ^1x1^2 + /3l2|V^l|V2 + /322|V'2r^2 + AV^i, 



1 72 

^ V^i(O) = ?/;i,o, ^'2(0)= ^2,0 

with X G M.^ and < 3. The parameters are the magnetic trap strength 7 > 0, the 
scattering lengths /3i2, /322 ^ and the external driven field constant A G M, called 
Rabi frequency. 

The energy and the mass associated to system (1) are the following 

m = I (^ivz^ip + ^iv^2p + ^ixp(iv^ip + i^2p) 

(2) + /3i2|^inV^2p + ^-f\^2\' + 2A3f?(^t^2)) {x.t)dx, 



M{t) = / (\ij,\' + \^2\'){x,t)dx, 



and are conserved along the flow of solutions to system (1). 

The above system models a mixture of Bose-Einstein condensates consisting of two differ- 
ent hyperfine states of Rubidium atoms confined in the same harmonic trap [3] . The first 
experiment concerning the binary Bose-Einstein condensate was performed at JILA with 
\F = 2,mf = 2) and |1, —1) spin states of ^^Rb [26]. By applying a weak magnetic (driven) 
field with Rabi frequency A, the two components were coupled in the overlap region. This 
coupling realizes a Josephson-type junction and gives rise to nonlinear oscillations in the 
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relative populations [37]. The effect of the Rabi term was examined in several papers 
[13, 27, 28, 33]. In [13, 28] it is shown that the system may be transformed into a system 
without the Rabi term given that all the interaction strengths are the same (the /3jj's are 
equal). Nistazakis et. al. [28] showed that due to the Rabi frequency it is possible to 
transfer the wave function of one component to the other. 

An interesting question in this framework regards the dynamics of a two-component Bose- 
Einstein condensate with mixed attractive/repulsive nonlinearities. Indeed, depending on 
the intra- and inter-species scattering lengths, the system will either reach a stable config- 
uration, or there will be occurrence of a singular dynamics. It is a standard procedure in 
physical experiments for single-component Bose-Einstein condensates, to use a Feshbach- 
resonance technique in order to tune the scattering length and to consequently stabilize the 
condensate. In mathematical language this is translated by introducing a time-oscillating 
coefficient in front of the nonlinearity in the Gross-Pitaevskii equation [1, 20, 31]. This 
stabilization method attracted a wide interest in physical literature and it is studied ex- 
tensively from the experimental, numerical and theoretical point of view. However, the use 
of Feshbach resonance would imply the strength of the applied magnetic field to oscillate 
at frequencies much higher than the characteristic frequencies of the system. To overcome 
this difficulty, Saito et al [32] proposed to use two hyperfine states with different scattering 
length. In such a system the Rabi oscillation between the states will cause an effective os- 
cillation of the scattering lengths and stabilize the system; this is performed by irradiating 
the system with a constant electromagnetic wave. Our purpose in this paper is to give a 
rigorous analysis for such systems. We investigate the dynamics of a two-component non- 
linear Schrdinger system with focusing and defocusing nonlinearities driven by an external 
field with Rabi frequency A. A considerable amount of mathematical literature is dedicated 
to the study of long time behavior of nonlinear Schrdinger systems (see references below), 
however the case of systems with a linear coupling term coming from the Rabi frequency 
is poorly understood. Numerical experiments performed in [17] suggest the fact that the 
Rabi frequency may affect the long time behavior of the system. Therefore here we address 
to the question whether the presence of an external driven field can prevent blow-up in 
the presence of focusing nonlinearities. More in particular we investigate the asymptotic 
behavior of the system when the Rabi frequency becomes large, and we rigorously study 
the limit as |A| — ?■ oo. We show convergence in Strichartz spaces of solutions to system (1) 
towards solutions to 

{1 7^ 
1 7^ 
idtU2 = --AU2 + —\x\'^U2 + <Jl\u2\'^U2 + (T2\ui\'^U2, 

with ai = {{3n + 2/3i2 + /322)/2, aa = {/3n + /322)/2, when |A| ^ oo. 

Moreover, the solution of the system (1) is proven to exist on a time interval strictly less 
than the maximal existence time of the limiting system. This implies in the case of global 
existence of the limiting solution, that the solution of the original system exists on finite 
time [0,T], with arbitrarily large T. Therefore, at least for large A we expect the original 
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nonlinear system to behave like the limiting system, which is nothing but a well-known 
nonlinear Schrodinger system. 

Here we give a short overview on the papers dealing with global existence/blow- up of the 
system. The following Schrodinger system without external potential and with power-type 
nonlinearities was extensively treated in the literature: 

idt^Pi = -^Av^i + (/^iiiV'ir^ + /3i2i^ir-V2r+')^i 

Fanelli and Montefusco [14] deal with this system with P12 < 0, and Pu = P22 = "l- 
Using the existence of ground-state solutions they establish a vector-valued Gagliardo- 
Nirenberg inequality, which also gives the mass threshold for global existence of solutions 
to the system in the mass critical case p = 2/N. Li, Wu and Lai [21] improve partially the 
arguments of Fanelli and Montefusco in the critical and super critical case p > 2/N. They 
obtain a new blow-up threshold in terms of the ground state. The same results for more 
general coefficients of the nonlinearities are established by Ma and Zhao [25] , by using the 
sharp vector- valued Gagliardo-Nirenberg inequality they derived previously in [24]. The 
planar case = 2 is treated by Ma and Schulze [23]. Lin and Wei [22] prove simultaneous 
blow-up of the two components by studying the associated self-similar solutions for p = 1. 
Song [34] considers the system with more general power-type nonlinearities and establishes 
the criterion for blow-up and global existence in the critical case. The blow-up profile is 
analyzed by Chen and Guo [11]. In the critical case p = 2/N they show that the solution 
blows up like the Dirac 6 distribution. Moreover, they carry out the result for the system 
with harmonic potential by using the Lens Transform [7, 35]. Note that in the case p = 1 
we recover (1). 

Among the papers dealing with nonlinear Schrodinger systems with an external potential 
we mention [12, 30, 38]. In [12] the authors consider (1) with all f3ij = —1 and derive a 
sharp criterion for the energy norm for blow-up / global existence using the existence of the 
ground state solution of the scalar equation without potential. Prytula et al. [30] study 
a system of M equations and prove a sufficient condition for blow-up depending on the 
coefficients Pij. The two-dimensional case N = 2 for a system (1) with harmonic potential 
was considered by Zhongxue and Zuhan [38]. Here the ground state solution of the system 
is related to the threshold condition for blow-up. 

There are only few papers considering the system with the Rabi frequency. Bao and Cai 
[4] treat the stationary case, and prove existence and uniqueness of ground states. In this 
framework they also perform asymptotics for |A| — t- 00 and obtain the limiting behavior 
of the ground states. The solitary wave solutions of a two-component system when all 
the physical parameters are spatially dependent is studied in [5]. Zhongxue and Zuhan 
[39] prove finite-time blow-up of solutions for initial data greater than a threshold when 
all nonlinearities are focusing. The results seems to be valid in the two-dimensional case 
only. The nonlinear Schrodinger system of two equations with an external driven field 
and with nonlinearities which may not be cubic is analyzed in [17]. The authors derive 
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sufficient conditions for blow-up and global existence. A semi-implicit formula shows the 
mass transport between the two components. 

In [10] Cazenave and Scialom study rigorously the asymptotic behavior of a nonlinear 
Schrodinger equation with a time-periodic coefficient in front of the nonlinearity. Fang 
and Han [15] prove the convergence in the case of an energy critical nonlinearity, thus the 
power of the nonlinearity is 4/(A^ — 2). 

This paper is organized as follows: in Section 2 we introduce the main notations and 
tools used throughout the paper. We furthermore study system (1), establishing the cases 
when there is global existence of solutions or possible occurrence of finite time blow-up, 
according to the choice of the coefficients /3jj's. Then in Section 3 we transform (1) in a 
suitable way, in order to study the system when the Rabi frequency A becomes large. We 
then write down the limiting system, and in Section 4 we prove the rigorous convergence 
result, namely Theorem 11, which is the main result of our paper. As the reader will see 
in Section 4 the proof of Theorem 11 follows the same steps as the paper by Cazenave and 
Scialom [10]. 



In what follows C will denote a generic constant greater than 1, which may possibly 
change from line to line. With 3? and S we denote the real and imaginary part of a 
complex number, respectively. By z* we denote the complex conjugate of z. The scalar 
product between two vectors v, v' will be denoted by {v, v'). 

Since we are dealing with two-component Schrodinger systems, we shall indicate with 
capital letters the two-dimensional vector fields describing the wave-function of a two- 
component quantum system. For example we shall write \E'* = (^/'i, ■?/'2); or \&q = {ipi^Q., ip2,o), 
^* = (^]L!^2) so on. Consequently, we also denote 



2. Notation and preliminary results 





where A, B[^] are the two matrices defined by 



(5) 



(4) 




The energy and mass may be written 



(6) E{t)= + ^\x\^\%^ + -^>*B[%^! + 2\^{^l)li)2)){x,t)dx. 
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We shall denote by Lp(M^), W^'P(R^), the usual Lebesgue and Sobolev spaces, respectively. 
We shall also make use of the mixed space-time Lebesgue (or Sobolev) spaces, so that for 
example L«(J; denote the space of those functions / having the following norm 

finite, 

1 

||/||l.(/;L^(m-)) := y^(^I^Jf{t,x)\'dxY^ dt 

We often shorten notation L^L^.^! x M^) = L^'(R^)) or LfL^. if there is no source of 
ambiguity. 

We are interested in studying system (1) in the energy space, which is defined by 

E(M^) := {/ e H\R^y, I ■ 1/ e L^{R^)}. 

It is straightforward to see that, for any \E' G S(]R^), the energy in (2) is finite. 
Let us consider the Hamiltonian with confining potential, H = — |A-|-^|xp, the associated 
propagator is So{t) := e"**"^. It is unitary on L^(M^) and from Melher's formula (see for 
example [8, 29]) one can see that it satisfies a dispersive estimate for short times, namely 

N 

(8) ||5'o(t)/|Uoc(i,iV) < |tr^||/|Ui(MiV), \t\ < 5, 
for some small 6 > depending on 7. 

Furthermore, in what follows we also need the following commutation formulas for H: 

(9) [V,H]=j^x, [x,H] = V. 

Even if (8) holds only for short times, we know the group So{t) enjoys the same Strichartz 
estimates as the free Schrodinger propagator without confining potential, but only locally 
in time, i.e. the constants appearing in the inequalities below depend on the length of the 
time interval. Here we state the results we need, for further details we address the reader 
to [8], Section 2. 

Definition 1. We say (g, r) is an admissible pair if 2 <r < (2 < r < 00 for N = 1, 

2 < r < 00 for N = 2), and 

1 iv /I r 



(10) 

^ ^ q 2 \2 r 

Proposition 2. Let {q,r), {q,r) be two arbitrary admissible pairs. Then, for any compact 
time interval, we have 

\\So{t)f\\LjLr{Ixm) <C(|/|,r)||/||i2(IRiV), 



/ Soit - s)F{s)ds\\L.Lr^i,^N) <C(|/|,r,f)||F| 
Jo 



The estimates in Proposition 2 will serve us to show local well-posedness of (1) in S(]R^). 
As we already said in the Introduction, in order to study the asymptotic behavior of (1), 
we first transform the system into a similar one, see (22), and then we perform the limit, 
obtaining (27). All those systems are quite similar and at a local level they can be studied 
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in the same way. For this reason here we write a more general local well-posedness result, 
which may be applied to all those systems involved in this paper. The result in Proposition 
below is quite standard; its proof consists in just adapting the method by Kato introduced 
in [18] for NLS equations, to the case of systems and for Hamiltonians with a confining 
potential. As it is already well explained in [9] (see Chapter 4 and also Remark 3.3.12), 
the known theory for nonlinear Schrodinger equations is easily extended to systems. In [9], 
[18] the presence of a confining potential is not considered, however this is only a minor 
modification for the local well-posedness framework, thanks to Strichartz estimates stated 
in Proposition 2. 

Proposition 3. Let us consider the following Cauchy problem 



where the unknown ^ is a complex vector field. Let M G C(C^; C^) he such that A/'(0) = 0, 
and M{^) = 7Vi(*) + Af{'^2), where Afi,Af2 € C(C2; C^) satisfy 

• iM(^) - ^2{^)\ < + i^ni^ - ^1, 

for all^,^ e C\ 

For any "^0 G S(]R^), there exist S = 5(||\E'o||s) > and a unique solution"^ G C([0, 5]; S(]R^)) 
to (11). Moreover we have 

(12) ||^||l°°([o,<5];E(r'V)) < 2C||\E'o||e(r^)- 

Furthermore, the solution can be extended to a maximal interval [0,Tmax), o,nd the blow- 
up alternative holds true, namely ifTmax < oo, then 



Finally, for any < T < T^ax o-nd any admissible pair {q,r), we have \E^, V\E', | ■ I^E' G 



(11) 



idt"^ = — 
^(0) = 



|A$ + ^|x|2*+AA(^), 



Irni ||V^(t)||L2(iKiV) = cx); 



L'?([0,T];L^(M^)). 



Proof. Let us define the space 



K := s.t. V*, I ■ I* G L^LlilO, 6] x R^)) n L^/^ LI{[0, 5] x M^)), 




where M, 6 > will be chosen later, endowed with the distance 



c/(^,^) := 11^ - ^ 



L-L2([0,<5]xM^) + 11^ - ^IIl«/^l4([o,5]xM^)- 
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By the hypotheses on A/'(\l/) we have 

IIM(^) -M(^)||i2(KiV) <C||^-^||i2(KiV) 

\\^^2i^)-^^2imL*/HR^) <cm\\l,^^N^ + ||^'|ii4(i,iv))||* - ^hn^^) 



V ] Af2m\\L^/s^m)<C\\'^\\%iw^)\\i V 



Let us now consider the operator 

G[^] := So{t)^o - i f So{t - s)Af{^){s)ds. 
Jo 

By using the commutation rules for the Hamiltonian H and Strichartz estimates, we have 



<c||vi/ob + c7|| I V ]^^^m\L}L^ + c\\\ V \Af2m\ 



8 

, H^TTV 4/3 



Analogously we obtain 



1 2 I 1 1 ,Tt 1 1 2 



By the Sobolev embedding and by using Holder's inequality in time in the 

previous expressions, we have 

/ 1 \ 

\\Gm\\L^EAmx^-) + \\ ( V jGmhrLUm.mn <C\\%h + G{6 + 6'^M')M 

||G[vI/i]-G[vl/2]||^s/.^4(p^,j,^,) <G{6 + 6'^M')d{^,^). 
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Now, we choose M, 5 such that 

^ II y^ollE(IR^) — 

so that we have 

l|G'[^']IU^([0,5];E(R-)) + ll ( V \G[mLrLiim^^-)^^ 



(13) 



d(G'[vl/],G[^'])<irf(vl>,^'). 



This imphes G is a contraction in therefore there exists a unique ^ & K such that 

= = 5o(t)*o - ? / Soit - s)M{^){s)ds. 

Jo 

Hence ^ G C([0, 5]; S(M^)) n ^^/^([O, 5]; ^^(M^)) is a solution to (11) in [0, 6]. From (13) 
we also see 

||^IU°°{[0,<5];E{I8'V)) < M = 2C||\&o||e(I8'V), 

which proves (12). Analogously, by Strichartz estimates we also have 
1 \ 

I ^IIl?lj([o,5]xm^) < C'll^olls + C(5^(i + M^)M < M. 

Furthermore, from the proof in the fixed point argument above we also infer that we may 
extend the solution as long as the L^— norm of the gradient of the solution remains bounded, 
hence the blow-up alternative holds true. This implies we can extend the solution to a 
maximal interval [0,Tmax), and moreover for any < T < T^ax, {<l,f') admissible pair, we 
have 

1 

^IIl«L5([0,T]x]R^) < oo. 

□ 

Corollary 4. The system (1) is locally well-posed in S(]R^). Furthermore the energy and 
the mass defined in (2) are conserved in [0,Tmax)- 

A natural question which arises at this point is to see whether the solution in Corollary 
above is global or there is possible occurrence of blow-up in finite time. 
For the case of a single cubic nonlinear Schrodinger equation the picture is complete. 

• If the nonlinearity is defocusing (i.e. its coefficient is positive) then, by using the 
conservation of energy and the consequent uniform a priori bound on the if ^— norm 
of the solution, global well-posedness holds true. 
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• For = 3, if the nonlinearity is focusing (i.e. its coefficient is negative), then there 
exist initial data for which the L^— norm of the gradient of the solution blows up 
in finite time. 

• In the case N = 2 and focusing nonlinearity, we have 

— if II'i/'oIIl^ < ||Q||l2, where Q is the unique positive radial solution to AQ + 

— Q = 0, then the solution exists globally [36]; 

— for initial data HV^oIIl^ > ||Q||l2 there is possible occurrence of blow-up. 
Unfortunately, in the case of system (1) the picture is not complete and we can give only 
a partial answer. 

Theorem 5 ([17]). Let N < 3 and set (3 = max{2(-/3n) + , 2(-/322) + , (-/3i2) + }- Then 
there exists a global-in-time solution to in the following cases: 



(1) /322, /3i2 > or < /3n^22 with /3n > 

(2) iV = 1 

(3) iV = 2 and M(0) < 2/(C2/3), tf mm{(3n, f322, M < 

(4) N = 3, ||V^(0)||2 < 2{E{0) + \X\M{0)), and M(0)(E(0) + |A|M(0)) < zf 
min{/3ii,/322,/3i2} < 

Proof. We show global existence for /3^2 < /3ii/322 with /^n > 0. It follows that 



is always positive, thus using the energy conservation we get uniform bounds on the S-norm 
of the solution 



In the next Theorem we show the cases in which there is possible occurrence of blow-up 
in finite time. We describe two situations, the first one in which the nonlinearity is negative 
definite (see below). In this case we apply a method introduced by Carles [6] for focusing 
NLS equations with a confining potential, by using a modified energy functional (for more 
discussions about this modified energy and its interpretation see [6]). In the second case we 
assume that at least one of the coefficients /3j/s is focusing and that the energy is negative 
enough so that conditions (ii) or (iii) below are fulfilled. In this case we apply the method 
by Glassey [16] using virial identities. 

Remark 6. Note that the following sufficient condition for blow-up is valid also in the 
supercritical case N = 3, whereas in [17] the authors consider only the mass-critical case 
N = 2. 

Theorem 7. Let G C([0, T^ax); ^(1^^)) be the solution to (1) as in Corollary 4, and let 
us define the virial potential 





3ft(V^*^2)(a;,t)rfx < E{t)+\\\M{t) = E(0)+|A|M(0). 



For all other cases see [17]. 



□ 
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Let us assume N >2 and one of the following conditions is satisfied 

(i) the nonlinearity is negative definite, i.e. we either have /3^2 ~ /3ii/322 < with 
< 0, or /3i2, /322 < 0, and we also assume that 



E(0) + ^M(0)<^/(0); 



(ii) min{/5ii,/322,/9i2} < 



(E(0) + |A|M(0))<^J(0); 



(iii) min{/3n,/322,/3i2} < 0; /'(O) < and 

{E{0) + |A|M(0)) < I'jO). 



iV + 2 ' ' ' V^Tn 

Then the solutions blows-up at a finite time, i.e. 3 < T* < oo, such that 

lim ||V^f(t)||L2 = oo. 

Proof. We first consider case (i). Analogously to [6] we introduce the following functional 

Ei(t) =cos^(7t) / (-|V^P + -^*S[^]^ + 2A3fJ(^/^i*^/;2))(t,a;)cia; 
Jj^N 2 2 

+ sin2(7t) / + ^sin(27t) / x-J{t,x)dx 

^ JRN 

where J{t, x) is the current density defined by 

J(x,t) = ^3(V^^V^/'l + V^2*V^^2)• 
Let us first notice that the condition E{0) + i^M(O) < ^/(O) implies Ei{0) < 0. By 
computing the time derivative of Ei we obtain 

= 7sin(27t)^^ / (t, A|7 sin(27t) / (2 sign(A)3f?(7/'*7/;2)-|*P) (t, 

By the assumptions on the coefficients /3j/s we have \E'i?[^]^ < for any 7^ 0, we thus 
infer 



(14) 1^1^^)-°' 



IT 

0,— 

. 27, 



On the other hand, we see that if exists in [0, then 

, = — / \x\ \ 
V27; 



E, (—] = ^ / \x\''\-^{x,t)\''dx > 0, 



ASYMPTOTIC BEHAVIOR OF NONLINEAR SCHRODINGER SYSTEMS WITH LINEAR COUPLINGll 

contradicting -E'i(O) < and (14). Thus, there exists T* < ^ such that 

hm ||V^'(t)||L2 = oo. 

For cases (ii) and (iii), we use the virial identities. By similar calculations as above we 
have 



I'{t) =2 X- J{x,t)dx 

I"{t)= I (2|V^P + iV^*S[^]^-272|xn^p)(t,x)d 



We write 



I"{t) = 2NE{t) + 2{2-N) I \Vm{x,t)\^dx-{2 + N)-i'^I{t)-AN\ I ^{iPli)2){t,x)dx. 

From the conservation of energy and mass we then obtain 

I"{t) = -(2 + N)-f^I{t) + 2N{E{0) + |A|Af(0)) + R{t), 

where 

R{t) = 2{2 - N) f \V^{t,x)\^dx-4:NX f ^{iPli)2){t,x)dx - 2N\\\M{t) <Q. 

The solution of the above differential equation is given by 

I{t) = cos(V2 + A^7t)/(0) + ^ ^ sin(V2 + N-ft)I'{0) 

V 2 + A^7 

2N 

+ 77^(^(0) + |A|M(0))(1 - cos(v/2TiV7t)) 



(2 + iV)72 

+ /TT^ / MV^ + N^it - s))R{s)d 
V2 + iV7 Jo 



< cos(V2 + A^7t)/(0) + ^=^^ sin(V2TiV7t)/'(0) 

V2 + A^7 

2 TV 

(E(0) + |A|M(0))(1 - cos(v/2TiV7t)). 



(2 + iV)7 



We then notice that, if (ii) or (iii) are satisfied, then J(t) will eventually become nega- 
tive, giving a contradiction. The unboundedness of the gradient then follows from the 
conservation of mass and the uncertainty inequality 

mmh < ■ i^(t)iuHiv^(t)iu. < ^)\\vm\\L^- 



□ 
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3. The transformed system 

As we showed in Theorem 5 the picture about global existence or possible blow-up of 
solutions to (1) is not complete. In order to get a better understanding, at least in the case 
when |A| becomes large, we are going to study the asymptotic behavior of solutions to (1) 
and the limiting system. In this Section we transform the system (1) into a similar one, 
which is more suitable to handle in order to study the asymptotic limit. Let us consider 
the linear part of (1), 

1 '7^ 

(15) i^t^ = --A^ + y|xp^ + A*. 

If A = (i.e. A = 0), then the two equations in (15) decouple and they evolve independently 
through the Hamiltonian H := — |A + and we have \E'(i:) = So(t)'^Q. In the case 

A 7^ the two equations are coupled, and thus 

where Sx{t) = So{t)nx{t) = e-''^Vtx{t), and 

O (+\ — ^-itA - f cos(At) -isin(At) 
ux[i).-e - 1^ _,sin(At) cos(At) 

It is straightforward to see that Qx(t) satisfies the following properties 

• Qxity = ^x{-t) for Vt G M 

• Qx{ti)^xit2) = Qx{ti+t2) for every ti,t2 e M. 

System (1) can now be written in the following integral formulation 

(16) ^(t) = Sxmo - I [ Sx{t - s)Bm^{s)ds. 

Jo 

To study the asymptotic limit when |A| — > oo, we want to cancel out the oscillations in the 
homogeneous part. We write 

(17) Qxi-mit) = Somo-^ [ Soit - s)Qxi-s)Bmnxis)Qxi-s)^{s)ds. 

Jo 



Now we define 



1 / 1 A 1/1 



and ^ ~ 75 1 ''^j^ ^ to be the matrix having vi,V2 as columns, so that Qx{'t)vi 
e~*'^*Wi, ilxit)v2 = e^^^V2-, and 



;i8) nx{t)V = V 
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By using Vi , f 2 we define 

so that ^(t) = e-'^^(f)i{t)vi + e'^^<p2{t)v2 and 

Vt{-t)<^{t) = (l)i{t)v^ + (I)2{t)v2. 

We want to find the system of NLS equations satisfied by $ = (0i,02)*- From formula 
(17) we can see that 



(19) 



$(t) = 5'o(t)$o - I [ So{t - s)V'nx{-s)B[^]nx{s)V<^{s)ds 
Jo 



where $0 '■= ^*^o- Now we notice that B[^] 
write 



^1 



and we may 



(20) 



where I_ 



^1 



^2 



1 
71 



(e-*^Vi(t)I + 6*^*02(^)1-) 



We want to compute the nonhnear potential matrix V^Qx{~s)B[$]Qx{s)V =: in 
formula (19). By using (18) and (20) we see it equals 



iXs 



V' (5|0i|2 + I_M_|02|' + 



h<P2 



2iXs 



I_50t02) V 



-iXs 



After some tedious calculations, and by rearranging a bit the terms we write 
(21) 



m =4 



(/Sll + 2/3i2 + /322) 



+ |(/3ll-/322) 



-2iXti 1* I „2iAt 



+ (/3n-2/3i2 + /322) 



-2iXt 



(l0iP + l02p) 



+ -(/3n-2/3i2 + /322) 



where we divided the autonomous part 13°° [$] from the non-autonomous one R'^ [$] . Con- 
sequently, the equation for the transformed variable $ becomes 



(22) 



1 -7^ 
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with B defined above in (21). Let us notice that in (22) the term i?°°[<l>]$ can be also 
written in the following compact form 



4 1 2(/3n + /322)0^02 (/3ii + 2/3i2 + /322)|02p ' ' 



1 / (/3ii + 2/3i2 + /322) 101' 2(/3n + /322)0i02 



thus we may rename to be the matrix in the expression above, i.e. 

/OON Roor^,! _ 1 (/3ll + 2/3l2 + /322)|0P 2(/3n + /322)0102 

^ ^ ^ ^ 4V 2(/3n + /322)0*i02 (/3ii + 2/3i2 + /322)|02p 

It is straightforward to check that the nonlinear terms B°°,R'^ satisfy 

(24) - B^[F2]F2\ < {\Fi\^ + \F2\^m " ^2!, 

(25) |V (^B^[F,]F, - B^[F2]F2^ \ < (|Fip + \F2\^)\V{F, - F^)] 

(26) \V\R^[F,]Fi)\ < iFinV^'Fil, = 0, 1 

Remark 8. Since the T.-norms o/\E' and of^^ are equal, all the results stated in Corollary 
4 and Theorems 7 and 5 hold true also for . 

Formally, in the limit as |A| goes to infinity, we expect the non- autonomous part to cancel 
out, because of the oscillating coefficients in front of the nonlinearities which average to 
zero. We thus obtain the following system of coupled nonlinear Schrodinger equations 

(27) idtU = --/\U + ^\x\^U + B^[U]U. 

and the corresponding limiting energy and mass are given by 
(28) 

m= [ (l-\VU\' + ^\x\'\U\' + lu*B^[U]u)ix,t)dx 
Jrn \z I I / 

|2 , Al + 2/3l2 + /322 ,4 ,1 |4n , /3ll+/322| ,21 \i\ I +\a 

I H ^ (Fi| + P2I )H ^ \u\\ \U2\ j{x,t)dx, 

M-{t) = / \uj{x,t)\'^dx, with j = 1,2. 

Remark 9. Let us consider formula (17) again, 

Qx{-t)^{t) = Somo-^ [ So{t-s)Qx{-s)BmQx{s)Qx{-s)^{s)ds. 

Jo 

If we define 

l>(t) :=fiA(-t)^(t), 
then we obtain the following system 

ft 



\vu\ 




' 1 r 


\vu\ 




7^, ,: 
2 



^t) = Somo-t f So{t - s)B2[^Ms)ds, 
Jo 
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but now the nonlinear potential matrix B2 [$] is given by the more complicated expression 



52[$]=^(/3ll+/522)|$| 



2 



1 

1 

1 - 

1 - 



(|01P-|02P) 



1 1 
-1 -1 



h92 



-2iXt 



+ e 



2i\t 



g(/3ll-/322) 
^(/3ll-/322) 



1 1 
-1 -1 

1 -1 
1 -1 



(|0ir+|02r) + ^(/3ll-/322)I010; 



(|01P+|02p) + T(/3ll-/322M02 



4 



1 



+ ^(/3n-2/3i2 + /322)<|e"^^"* 



1 1 

-1 -1 



>W2 



>lV2 



Furthermore, if we only consider the autonomous part, we see that even in the asymptotic 
limit we obtain an expression for the nonlinear potential matrix which is rather compli- 
cated. 

On the other hand, if all the inter- and intra-species coefficients equal, i.e. = P12 = 
(322 = P, then the expression (and the analysis) simplifies considerably. Numerical studies 
using this transformation in the case of equal coefficients are performed in [13, 28]. Never- 
theless our aim in this paper is to consider a general choice for the coefficients Pij 's, that 
is why we choose to transform the system in terms of ^. 



4. ASYMPTOTICS FOR A OO 

In this Section we prove the rigorous convergence of solutions to system (22), towards 
solutions to system (27), when the parameter |A| goes to infinity. As the reader can see it 
is heavily inspired by the paper by Cazenave, Scialom [10]; here a system of NLS equations 
is considered, and we also have to deal with the presence of a confining potential, but the 
strategy of proof follows the same steps as [10]. 

The main result of this Section is Theorem 11, but before stating it we first analyze system 
(27). By Proposition 3 and Theorem 5 we have the following 

Proposition 10. Let Uq G S(M^). Then there exist 6 = 6{\\Uo\\j:(m)) > and U G 
C([0, 5]; S(]R^)) unique solution to system (27). Moreover, 

(29) l|f^llL°°([0,(5];S(M^)) < 2C||f/o||s(R'V)- 

Furthermore, the solution U can be extended to a maximal interval [0,Smax), we have 
U G C([0, ^max); S(]R^)) and the usual blow-up alternative holds, i.e. if S max < 00, then 

\\VU(t)\\L2 ^00, ast /"Smax- 

Moreover, Smax = 00, thus the solution U to (27) is global-in-time in the following cases: 

• + /322 > and /3n + 213^ + P22>0 

• + /322 < and f + /322I < /3i2 
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For any < T < S^ax, admissible pair, we have 

U,VU, \ ■ \U e L'^{[0,T]; U{R^)). 

Finally, the three quantities defined in (28), namely the total energy and the mass of each 
component, are conserved in [0,Smax)- 

Now we can state our main Theorem. 

Theorem 11. Let e S(M^). For any X E M., we denote by the unique maximal 
solution to (22). Let U be the solution to (27), with initial data f/(0) = $0; in [OjS'max)? 
where < Smax < oo, as in Proposition 10. 

• For any < T < Smax, the solution exists in [0,T] provided \\\ is sufficiently 
large. 

• For any < T < Smax, admissible pair, we have 

(30) 11$^ - U^.^r + ||V($^ - U)\\^Lr + III ■ |($" - f/)|U?LS ^ 0, 

as \X\ — > oo, where the LfL^.— norms are taken in the space-time slab [0,T] x M^. 
In particular, convergence holds in C([0,T]; E(]R^)). 

The idea of the proof can be described as follows; let us consider the equation satisfied 
by the difference — U. By using the integral formulation we may write 

$^(t) - U{t) =-i ! So{t- s) - B°°[U]U] {s)ds 

Jo ^ 

-I I So{t- s)R^[<^^]^^{s)ds=: -ih-ih, 



see (21), (23) for the definitions of 

The oscillating coefficients in i?'^ converge weakly to zero, and since they appear inside the 
time integral in the Duhamel's formula, then I2 converges (strongly) to zero. For the Ii 
part, on the other hand, we can use a Lipschitz estimate for the nonlinearity B^[F]F and 
close the convergence with a continuity argument. 

To prove Theorem 11 above we proceed in two steps: first we prove that, as long as we 
have uniform bounds on on a space-time slab, we obtain the convergence. Then we 
prove that we indeed have such bounds on in [0, T] x M^, with < T < Smax, provided 
|A| is sufficiently large. 

We start by giving a technical Lemma which will be used to prove the convergence of I2 
to zero. 

Lemma 12. Let {q,r) be an admissible pair, <T < 00, and letfe L^' {[0,T]; L^' (R^)). 
Then, for any admissible pair {q,r), we have 

(31) II / So{t - s)e'^^f{s)ds\\L9(io,T];LriM.N)) -> 0, as \X\ 00. 

Jo 
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Proof. First of all, let us notice that by Strichartz estimates we have 

So{t - s)e'^'f{s)ds\\L,([o^T];L-iRN)) ^ 11/11^9' 
* 

Then, it is straightforward to see that by using a density argument it suffices to prove (31) 
for all / G C^([0, T]; iS(M^)). Let us consider the integral in (31), by integration by parts 
we obtain 

e^^^Soit - s)f{s)ds = ^ (e^^V(t) - 5o(t)/(0)) - ^ So{t - s) [dj{s) - Hf{s)] ds. 
Thus again by Strichartz estimates we get 

II ^* Soit - .)e^^7(^)rf^llL?L5 < ^ (||/||l?ls + ||/(0)|U^ + \\dj - HfW^,^^ , 

and the right hand side goes to zero as |A| — ?■ oo. □ 

Remark 13. In (31) we proved the convergence for the term with e^'^^ , hut this clearly 
holds for any other periodic function which has zero average, in particular it also holds for 
all the coefficients appearing in front of the nonlinearities in R^, see (21). 

The next Proposition shows that, as long as we have uniform bounds on the L^Sj,— norm 
of on a space-time slab, then we can prove the convergence of $^ towards U in that 
space-time slab. 

Proposition 14. Let $o G S(]R^) and let denote the maximal solution of (22). Let U 

he the maximal solution of (27), defined on [0, Smax)- Le-t < ^ < Smax (^nd assume that 
$^ exists on [0, P\ x and that 

(32) limsup ||$^||Loo([o,<?];E(MiV)) < oo. 

I A|— >oo 



Then we have 



1 



(33) lim II V (<f^-f/)|U,([o,£];L'-(MiV)) = 

|A|— >oo 

for any admissihle pairs {q,r). In particular U m L°°([0, £]; S(R^)). 

Proof. By hypothesis (32), we can take L large enough such that 

sup ||'^>^||l°°([0/];S(M^)) < OO. 
|A|>L 

By using Strichartz estimates similarly to the proof of Proposition 3 we see that from the 
bound above we also infer 

(34) sup sup II I V I $^||l?lj([o,£]xr'V) < oo. 
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where the second supremum is taken over all admissible pairs (g, r). If is a solution to 
(22) and U to (27), then we have 

$\t) - U{t) So{t- s) - fi°°[f/]f/) {s)ds 

-I [ So{t - s)R^[^^]^\s)ds =: -ih - ih, 
Jo 

where B°°, i?'^ are defined in (21), (23). We first consider Ji; by using Strichartz estimates, 
inequality (24), Holder's inequality and then (32), we have 

I|A||l<7([o,^];L'-(M^)) <C'(ll*^llL-([0,£];L-4(M^f)) + \W\\l'^{[0/];LHK'^)))\\^^ ~ 



<C||<l>^-t/||, 



'L5^([0/];L4(IRiV)) 

'LS^([0,£];L4(K]V))' 

For I2 we use Lemma 12, indeed by Sobolev embedding we have 



and consequently 

' So{t-s)R^[<^']<^\s)ds\\L.LraoA> 

converges to zero as |A| — > 0, for any admissible pair {q,r). By choosing (g, r) = we 
have 

\\^^ -U\\,s/N,^,,^,, ^j^, < e + C\\^^ -U\\ s 

II IlL/ L4([o,^]xK^') _ II "Lf^L4([0/]xR^) 

A standard continuity argument gives us 

11$^ - U\\^8/N^,^^^^^^^^^^ < Ce, 

which proves the convergence ^ U in L^^^L^([0,£] x M^). The same holds for any 
admissible pair (g,r), since 



where in the last inequality we used Holder in time. Consequently — t/ in LjL^.([0, i] x 
M^) for any admissible pair. We now prove the convergence of V$'^ and | ■ |$'^ in the same 
spaces. Again we consider the difference 

V($^ -U)+ x($^ - U) =: -ih - ih - ih, 
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where, by using the commutator relations of the Hamiltonian H with x and V, see (9), we 
obtain 

■= So{t - s)V - B°°[U]U^ is)ds, 

/2 So{t - s)x (^5°°[<l>^]$^ - B^[U]U^ {s)ds, 

/3 := f So{t - s) (V(i?^[$^]$^) + xi?^[<l>^]$^) {s)ds. 
Jo 

By using Strichartz estimates, inequahty (25) and Sobolev embedding, we have 

Wl^hrmoA.m) < C {W^'Wl^m^ + WWl^Hi) l|V($^ - U)\\^^^^j 
Analogously, for I2 we have 

wiArLtmix^n ^ ^ (ii^'iil-Hi + mwi^H,) III ■ - 

Finally, for I3 we use inequality (26), Sobolev embedding and (34) to estimate 

||v(i?^[$^]$^)+Mi?^[$1$^||^^^^^^3 < mrLrm {W^'^XrLi + HI " i^'iLr^O < 

Hence we can apply Lemma 12 and infer that 

11-^3 ||l?l-([o,£]xr'V) 0,as |A| 00, 
for any {q, r) admissible pair. By resuming, we have obtained 

||V(<l>^-f/)||^s/.^, + ||M(<l>^-f/)||^s/.^, < e+C (||V(<I>^ - U)\\^s/.^, + III ■ |(<l>^ - U)\\^s/.^,) . 
As before, a standard continuity argument gives us 

||V($^ - f/)||^B/.^, + III ■ \{^^ - U)\\^s^.^, < Ce, 

which shows the convergence in L^^^ as |A| — > 00. The convergence in LjL^([0,£] x M^) 
for any (g, r) admissible pair then follows as before, by using Strichartz estimates. □ 

The Lemma above states that, as long as the L°°([0, i]; S(M^))— norm of the family {$'^} 
of solutions for the Cauchy problem related to (22) is definitely bounded (as |A| — 00) in a 
time interval [0,i], then the family converges in to U solution to (27) with the same initial 
datum, in every Strichartz space. 

It thus remains to prove that given any time interval [0,T] with < T < Smax, the 
L°°([0, T]; E(M^))— norm of is indeed uniformly bounded, if we take |A| sufficiently 
large. 

For this purpose, let then T be a positive time strictly less than the maximal existence 
time for the solution U of (27), i.e. < T < S^ax, and let us fix the constant 

M := 2C||f/||ioo([o,T];S(R^))5 
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where C is the constant in (29). Let 5 = 6{M) be the constant in Proposition 3. Then 
exists in [0, 6] for all A and furthermore 

sup ||$'*'||l°°((0,(5);E(R'V)) < 2C||$o||s, 

Asm 

by (12). Now, let < £ < T be such that exists in [0,i], and that we have 

limSUp ||$^||L°°((0,£);S(IR^^)) < OO. 
|A|— s>oo 

Notice that by the inequality above we see that we can always choose i = 6. Thus, by 
Proposition 14 we have 

lim (11$^ - ?7||l9((0,£);L'-(R'V)) + ||V($^ - f/)||L9((0,£);L'-(R'V)) + III ■ |($^ - f/)||L9((0,£);L'-(R'V))) = 
I A|^oo 

for all admissible pairs (g, r). In particular we have 

hm \\^\i)-umm-) = ^- 

|A|— >oo 



This implies that 



sup ||$^(£)||2(MiV) <M, 
|A|>A 



for some A > sufficiently large. 

We can thus repeat the same argument, starting at time t = i. Consequently we have the 
solution exists in [0, 5 + £] and we have a uniform bound for the S— norm of 

sup ||$^||l°°((0,5+^);E(R'V)) < 2 CM. 
|A|>A 

Thus we can apply again Lemma 14 in the time interval [i, 5+i] and obtain the convergence 
of to U. Again, because of the convergence we have 

lim \\^\6 + £) - U{6 + i)\\^ = 0, 

|A|— >oo 



and in particular 



sup ||$\5 + £)||E(RiV) < M. 
|A|>A' 



Thus we repeat this argument to prove the result in the whole time interval [0,T]. 
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Appendix A. Simple example 

To illustrate the application of the presented theory we give a very simple example. Let 
us consider the case of one focusing and one defocusing nonlinearity, where the interspecific 
scattering length is zero. Thus we consider the following system in = 2: 

(35) idt"^ = — A^ + — + + 

with A as defined in (4), and 

We choose initial data such that the condition (3) in Theorem 5 is not satisfied, thus the 
initial mass M(0) is not too small, and none of the conditions (i)-(iii) of Theorem 7 is 
satisfied. In this way, we are actually not able to say if the solution may blow-up or exist 
globally. 

Remark 15. IfX = the system is decoupled and we have two scalar nonlinear Schrddinger 
equations, for which we know that ip2 exists globally, and ipi may blow up in finite time for 
initial mass ||'i/'i,o|U2(K2) > ||Q||l2(m^) with Q being the unique positive radial solution to 
AQ + Q^ -Q = 0. 

Performing the asymptotics for |A| — )■ oo we observe that \I/ converges to U solution to 
the linear Schrodinger system: 

1 '7^ 

2 2 ' ' 

Applying Theorem 11 we can say that the solution ^' of (35) exists in [0, T], for any T < oo, 
provided that A is taken sufficiently large. 
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